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Abstract 

We consider double-scaling limits of multicut solutions of certain one matrix models that 
are related to Calabi-Yau singularities of type A and the respective topological B model via 
the Dijkgraaf-Vafa correspondence. These double-scaling limits naturally lead to a bosonic 
string with c < 1. We argue that this non-critical string is given by the topologically twisted 
non-critical superstring background which provides the dual description of the double-scaled 
little string theory at the Calabi-Yau singularity. The algorithms developed recently to solve 
a generic multicut matrix model by means of the loop equations allow to show that the scaling 
of the higher genus terms in the matrix model free energy matches the expected behaviour 
in the topological B-model. This result applies to a generic matrix model singularity and the 
relative double-scaling limit. We use these techniques to explicitly evaluate the free energy 
at genus one and genus two. 



1 Introduction 



In [1], the large N limit of a class of J\f = 1 supersymmetric U (N) gauge theories was studied. 
The theories are in partially confining phase where an abelian subgroup G of the gauge 
group remains unconfined. The large N spectrum contains the usual weakly interacting 
glueballs, which are neutral under G, and also baryonic states which are electrically and 
magnetically charged with respect to G, whose mass grows like N. The models studied 
include the /3-deformation of M = 4 Super Yang-Mills and M = 1 SYM coupled to a single 
adjoint chiral superfield with a polynomial superpotential. At some isolated points in the 
parameter/moduli space of the models, these baryons can become massless, and this causes 
the 1/N expansion to break down. However, it is possible to define a double-scaling limit in 
which N goes to infinity and the mass M B of these states is kept constant. The crucial feature 
of this double-scaling limit is that there is a sector of the Hilbert space of the theory which 
decouples from the rest and has finite interactions which are weighted by the double-scaling 
parameter 1/N e ff ~ VT/Mb, where T is the tension of the confining string. Furthermore, 
it was proposed in [1] that the dynamics of this emergent sector has a dual description given 
in terms of a non-critical superstring of the type introduced in [2]. This dual formulation 
has the virtue that the worldsheet theory is exactly solvable and that the background is free 
from Ramond-Ramond fluxes. 

The exact vacuum structure and F-terms of the M = 1 models with an adjoint chiral 
field and a polynomial superpotential can be analyzed by means of the Dijkgraaf-Vafa matrix 
model correspondence [3-5]. Indeed the proposal of [1] is based on a careful analysis of these 
F-terms. The breakdown of the 1/N expansion corresponds to a singularity of the matrix 
model spectral curve and therefore of the dual Calabi-Yau. The baryon states that become 
massless correspond to Z)3-branes wrapping shrinking 3-cycles in the Calabi-Yau. 

In [6], this analysis was extended to a more general class of singularities. Again, it was 
found that at these particular points in the moduli space certain states become massless 
and that in a suitable double-scaling limit, where the mass of these states is kept fixed, 
a particular sector of the theory emerges with interactions governed by the double-scaling 
parameter. There are two novel features in these models. First of all, contrary to the cases 
considered in [1], there is no supersymmetry enhancement in the double-scaling limit. This is 
signalled by the fact that the glueball superpotential does not vanish in the interacting sector. 
In fact, this is also one of the reasons why the dual string background is not determined 
explicitly. Secondly, in [6], some or all of the states that become massless are neutral under 
the abelian subgroup G of the U(N) theory which remains unconfined. As a consequence, 
the presence of these extra massless states may not affect the coupling constants of G but 
is captured by the higher genus terms of the matrix model free energy as in [7]. These 
terms control certain F-term interactions of the glueball fields with the graviphoton and 
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gravitational backgrounds [8,9]. 



Another important feature that emerges from the analysis of [6] is that these large N 
double-scaling limits correspond to double-scaling limits of the Dijkgraaf-Vafa matrix model 
of the same kind that was considered in the "old matrix model" era to study c < 1 systems 
coupled to two-dimensional gravity [10]. In particular, in [6] it was shown that the double- 
scaling limits have a well-defined genus expansion in the sense that the genus g free energy of 
the matrix model F g scales like A 2 ~ 29 ~ M^~ 2g [6]. On the other hand, the singularities and 
double-scaling limits considered in [1,6] generally fall into different universality classes from 
the ones usually considered in the old matrix model. It is natural to ask what is the bosonic 
non-critical string that corresponds to these matrix model double-scaling limits and what 
is the relation between the bosonic non-critical string and the non-critical superstrings that 
enter in the dual description of the models considered in [1]. The answer to the first question 
is provided by the Dijkgraaf-Vafa correspondence [3,11,12] that states that a generic one 
matrix model with superpotential is mapped to the topological B model on a non- 

compact Calabi-Yau of the form 

uv + y 2 + W'(x) 2 + deformations = (1.1) 

In taking a double-scaling limit, we tune the parameters of the superpotential and the 
deformation polynomial so that we are in the neighbourhood of a particular singularity 
of the above family of Calabi-Yaus. For instance in [1] we are close to an A n _x singularity 

uv + y 2 = x n - fi , (1.2) 

whereas for the (2, 2p + 1) bosonic minimal model coupled to 2d gravity we would have 

uv + y 2 + x(x — ei) 2 ...(# — e p ) 2 = . (1.3) 

Therefore, we conclude that the bosonic non-critical string corresponding to the matrix 
model double-scaling limit of [1] is the topological B model at an A n ^i singularity. The case 
n = 2 corresponds to the conifold singularity. 



A check that this is consistent is provided by the fact that the scaling of the matrix model 
free energy F g ~ A 2 ~ 2a matches exactly the scaling of the topological B model free energy 

/ r \2-2g 

F toPtg ~^jn) g>\ (1.4) 

where Q is the holomorphic 3-form on the Calabi-Yau [9]. In fact, the double-scaling pa- 
rameter A corresponds precisely to the holomorphic volume of the 3-cycles that vanish at 
the singularity. This is in turn proportional to M B , the mass of the baryonic states, which 
come from _D3-branes wrapping the shrinking supersymmetric 3-cycles. Furthermore, the 
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Figure 1: The bosonic non-critical string defined by the matrix model double-scaling limit 
at an v4 n _i singularity corresponds to the A-twist of the above SL(2)/U(l) x LG worldsheet 
theory. 



fact that in the double-scaling limit F g ~ A 2-29 is a general result, it does not depend on the 
particular class of singularities one considers. It was derived in [6] by means of the recent 
algorithms to solve matrix models based on loop equations [13,14]. These allow to consider 
more general classes of singularities than were previously accessible via "old matrix model" 
techniques. 

The non-critical superstring backgrounds that appear as dual to the large N double- 
scaling limits studied in [1] are of the form 

R 3 ' 1 x (SL(2) k /U(l) x LG(W = X n )) /Z n , k = —— , (1.5) 

n + 2 

where LG(W) denotes a Landau- Ginzburg theory with superpotential W. They were initially 
introduced in [15] as holographic duals to the 4d double-scaled Little String Theory (DSLST) 
at a CY singularity of type A n _i, generalizing the proposal of [16] and previous work [17,18]. 
The non-trivial part of the above background has central charge 

k+2 n-2 . . 

c = c si + c LG = — 1 = 3 , (1.6) 

k n 
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and it corresponds to the geometry [19] 

fxz' k + uv + y 2 + x n = , 



(1.7) 



where z,u,v,y,x are homogeneous coordinates. This is equivalent to (jl.2|) . 

We argued, using the Dijkgraaf-Vafa correspondence, that the matrix model double- 
scaling limits considered in [1] are equivalent to the topological B model at an A n _\ singu- 
larity (jl.2J) . On the other hand, the matrix model captures the F-terms or topological terms 
of the 4d DSLST [1] which are given by the topological sector of the SL(2)/U(1) x LG back- 
ground (|1.5|) [15]. Therefore, we expect the non-critical string defined by the matrix model 
double-scaling limits to be associated to a topologically twisted SL(2)k/U(l) x LG(X n ) 
background. 

This proposal fits nicely with certain known results about the topological twist of the 
above background in the conifold case, n = 2, where the LG model is trivial. In fact 
in [20], Ghoshal and Vafa argued that the A-twisted N = 2 SL(2)/U(1) supersymmetric 
coset describes the topological B model on a deformed conifold. In [21], Mukhi and Vafa had 
previously shown that the above A-twisted coset at level 1 is equivalent to the c = 1 non- 
critical bosonic string compactified on a circle at self-dual radius. The open and closed sides 
of this map were recently analyzed in [22]. Therefore, as a direct generalization of the conifold 
case, we expect that the non-critical bosonic string defined by the double-scaling limit of [1] 
at an A n -\ singularity should correspond to the A-twist of the above SL(2)/U(1) x LG 
theory. In particular, for n = 2, the matrix model double-scaling limit should be equivalent 
to the c = 1 string. This fact can be checked directly on the matrix model side. Indeed, in 
the limit, the matrix model spectral curve becomes equivalent to that of a Gaussian model [6] 
which is equivalent to the c = 1 non-critical string [3] . This particular singularity is obtained 
from a 2-cut solution with a cubic superpotential in the limit where the two cuts touch each 
other. The fact that this singular limit should be related to the c = 1 non-critical string was 
also observed in [23]. 

The relation between strings on non-compact Calabi-Yaus and non-critical superstring 
brackgrounds [15,19] involving the N = 2 Kazama-Suzuki SL(2)/U(1) model or its mirror, 
N = 2 Liouville theory [15,24,25], has been studied by several authors (see [26-28] and 
references therein). 

Furthermore, the relation between the topological sector of six- dimensional DSLSTs de- 
fined at K3 singularities, the dual topologically twisted non-critical string backgrounds which 
generalize ()1.5|) . and certain non-critical bosonic strings, the (l,n) minimal bosonic strings 
has been recently studied in [29-31] (see also [32-34] for related matters). 

In the paper, we are going to use the matrix model double-scaling limit to study the 
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relative topological B model and non-critical bosonic string. In section 2, we will review the 
matrix models studied in [1,6] and the respective double-scaling limits. In section 3, we will 
review the proof that the genus g matrix model free energy F g goes like A 2 ~ 2g as shown in [6]. 
As we said, this argument applies to a general matrix model double-scaling limit and shows 
that the scaling of F g is consistent with the expected behaviour of the topological B model 
free energy F g j op ~ (/ 2§ [9]. In section 4, we will evaluate the genus one free energy F\ 
at the A m _i singularities considered in [1]. This gives information on the states that become 
massless at the singularity. In section 5, we compute the genus two free energy relevant to 
the matrix models considered in [1]. The result shows concretely how the double-scaling 
limit of F 2 depends on the details of the near-critical spectral curve. In the conifold case, 
the near-critical curve is a Riemann sphere and the general expression simplifies drastically 
and matches the well-known result. 



2 The double-scaling limit 



In this section, we will review the matrix model singularities and relative double-scaling 
limits studied in [1,6]. Consider an M = 1 U(N) theory with a chiral adjoint field $ and 
superpotential W($). The classical vacua of the theory are determined by the stationary 
points of W($) 

e+i 

9i. 



W($) = NTr 



N 



^ i 



(2.1) 



The overall factor N ensures that the superpotential scales appropriately in the 't Hooft 
limit. For generic values of the couplings, we find t stationary points at the zeroes of 

i 

W\x) = Ne l[(x - ai ) , e = g e+1 . (2.2) 
i=\ 

The classical vacua correspond to configurations where each of the N eigenvalues of <3> takes 
one of the £ values, {a^}, for % — 1, . . . ,£. Thus vacua are related to partitions of N where 
Ni > eigenvalues take the value Oj with Ni + N 2 + . . . Ng = N. Provided iVj > 2 for all i, 
the classical low-energy gauge group in such a vacuum is 

Gd = n u ( N i) ~ n x su ^ ■ ( 2 - 3 ) 

Strong-coupling dynamics will produce non-zero gluino condensates in each non-abelian 
factor of G d . If we define as W ai the chiral field strength of the SU(Ni) vector multi- 
plet in the low-energy theory, we can define a corresponding low-energy glueball superfield 
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Si = — (l/327i 2 )(Tr N .(W a iW m )) in each factor. Non-perturbative effects generate a super- 
potential of the form [35-37] 

e i 
. . . , S t ) = N ASj Iog(A}/5,-) + Sj) + 2m £ , (2.4) 

3=1 3=1 

where the bj are integers defined modulo Nj that label inequivalent supersymmetric vacua. 

Dijkgraaf and Vafa argued that the exact superpotential of the theory can be determined 
by considering a matrix model with potential W(&) [3,4] 

/oo 
rf$ exp {-g- 1 Tr W($)) = exp ^ F g gf~ 2 (2.5) 

9=0 

where $ is an N x N matrix in the limit N — > oo. The integral has to be understood as a 
saddle-point expansion around a critical point where Ni of the eigenvalues sit in the critical 
point aj. Note that N is not related to the N from the field theory. The glueball superfields 
are identified with the quantities 

t 

Si = g a Ni, S = J2^ = 9sN (2.6) 

i=i 

in the matrix model and the exact superpotential is 

^ BF ^ 

3=1 3 3=1 

where Fq is the genus zero free energy of the matrix model in the planar limit. 
The central object in matrix model theory is the resolvent 

w(x) = i'&^- s . (2.8) 
N x - $ 

At leading order in the 1/N expansion, u>(x) is valued on the spectral curve E, a hyper-elliptic 
Riemann surface 

y 2 = ^ 2 (W\xY + f^(x)) . (2.9) 
The numerical prefactor is chosen for convenience. In terms of this curve 

u{x) = W'(x) - Ney(x) . (2.10) 
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In (12. 9|) . fi^i(x) is a polynomial of order I — 1 whose £ coefficients are moduli that are 
determined by the S{. In general, the spectral curve can be viewed as a double-cover of the 
complex plane connected by I cuts. For the saddle-point of interest only s of the cuts may 
be opened and so only s of the moduli fe-i{x) can vary. Consequently y(x) has 2s branch 
points and I — s zeros: 1 

E: y 2 = Z m {x) 2 a 2s {x) (2.11) 

where I = m + s and 

m 2s 

Z m (x) = - Zj) , (J 2s {x) = - <jj) . (2.12) 

3=1 3=1 

The remaining moduli are related to the s parameters {Si} by (|2.6|) 

S i =g s N i = Ne I ydx , (2.13) 
where the cycle Aj encircles the cut which opens out around the critical point a« of W(x). 



Experience with the old matrix model teaches us that double-scaling limits can exist 
when the parameters in the potential are varied in such a way that combinations of branch 
and double points come together. In the neighbourhood of such a critical point, 2 

y 2 — > CZ m {x) 2 B n {x) , (2.14) 

where Zj, bi — > Xo, which we can take, without loss of generality, to be xo = 0. The double- 
scaling limit involves first taking a —>■ 

x = ax , Zi = azi , bj = abj (2-15) 

while keeping tilded quantities fixed. In the limit, we can define the near-critical curve XL: 3 

£- : yi = Z m (x) 2 B n (x) . (2.16) 

It was shown in [6], generalizing a result of [1], that in the limit a — > 0, in its sense as a 
complex manifold, the curve S factorizes as X_ U £+. The complement to the near-critical 
curve is of the form 

S + : y\ = x 2m+n F 2s . n (x) . (2.17) 
where F 2s - n (x) is regular at a = 0. 

1 Occasionally, for clarity, we indicate the order of a polynomial by a subscript. 

2 We have chosen for convenience to take all the double zeros {zj} into the critical region. 

3 For polynomials, we use the notation f(x) — Yiii^ ~ fi)> where f(x) = Yli( x ~ fi) ■> x = ax anc ^ fi = a fi- 
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It is important to stress that the above singularities are obtained on shell [1,6]. The 
family of spectral curves ()2.14|) corresponds to vacua of a given field theory. As such, the 
family satisfies the F-term equations coming from the exact superpotential ()2.7|) relative to 
the model and the choice of semiclassical vacuum. The solution to the problem of engineering 
these singularities on shell, namely the problem of finding a field theory model and tree level 
superpotential whose spectral curve exhibits the desired singularity in its moduli space, is 
explained in detail in [6]. The case where there are no double zeroes, m = 0, has been 
studied in [1,38,39]. The tree-level superpotential can be taken to be 

W($) = Ne Tr N [$ n+1 - W$] , (2.18) 

and the relative on-shell spectral curve is 

y 2 = ( x n - uf - Ul . (2.19) 

At each of the critical values U = ±U C , n branch points collide and the curve has an A„_i 
singularity. For instance, as U ^U c 

y 2 ^x n -(U-U c ) . 



In the a — > limit, it was shown in [6] that the genus g free energy gets a dominant 
contribution from £_ of the form 

F g ~ (AT a ( m +"/ 2+1 )) 2 - 29 . (2.20) 

Note that in this equation iV is the one from the field theory and not the matrix model N. 
This motivates us to define the double-scaling limit [1,6] 

a ^ , N ^ oo , A = Na m+n/2+1 = const . (2.21) 

Moreover, the most singular terms in a in ()2.20|) depend only on the near-critical curve ()2.16|) 
in a universal way. 



Observe that Eq. ()2.20|) matches the expected behaviour of the topological B model free 
energy at the singularity [9]. In fact, as can be seen from (|2.14j) and ()2.15|) 

A ~ N J ydx . (2.22) 

More precisely, the double-scaling parameter is proportional to the period of the one-form 
y dx on one of the cycles that vanish at the singularity. Moreover, this one-form corresponds 
to the reduction of the holomorphic 3-form Q on the underlying Calabi-Yau geometry 

uv + y 2 = W'(x) 2 + f(x) (2.23) 
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n = dudvdx (2 24) 

y/uv - W'(x) 2 - f(x) 

This comes from the fact that 3-cycles in the Calabi-Yau correspond to two-spheres fibered 
over the complex x plane. In particular 

tt ~ J ydx (2.25) 

where Q is integrated on a vanishing 3-cycle in the Calabi-Yau that reduces to one of the 
vanishing one-cycles in the matrix mode spectral curve. Putting everything together, we 
find that 

( f \ 2 ~ 29 f r \ 2-29 

F g ~A 2 ~ 2g ~[ ydx) (2.26) 



which is precisely the behaviour we expect for the free energy of the topological B model 
on the Calabi-Yau [9], in agreement with the Dijkgraaf-Vafa correspondence. In section El 
we will review the proof of Eq. (J2.2(Jj) [6] using the algorithms of [13, 14] based on the matrix 
model loop equations and we will see relation ()2.22|) arise naturally. It is worth stressing that 
this result is general, in the sense that is does not depend on the specific kind of singularities 
one considers. In this particular respect, the methods used are more powerful than "old 
matrix model techniques" , where one is usually limited to considering one-cut matrix model 
solutions. 



3 The Matrix Model double-scaling limit: higher genus 
terms 



In this section, we will be concerned with the behaviour of the higher genus terms of the 
matrix model free energy in the limit a — > 0. The most powerful methods for calculating 
the F g involve orthogonal polynomials (see the reviews [10]); however, these techniques 
have only been successfully applied to the case when the near critical curve has at most 
two branch points (but any number of zeros). The only known way to calculate the F g in 
general involves analysing the loop equations and in particular using the algorithms recently 
developed in [13, 14]. In the following, we will review these algorithms and the proof based 
on them that 

F g ~ A 2 " 29 (3.1) 

which was given in [6]. As we said above, this is the behaviour we expect given the Dijkgraaf- 
Vafa correspondence between the matrix model and the topological B model. 
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3.1 The loop equations 



The p-loop correlator, or p-point loop function, is defined as 




X — $ / conn 



P > 2 . 



(3.2) 



It has the following genus expansion 



oo 1 
9=0 N 9 




(3.3) 



In [13], Eynard found a solution to the matrix model loop equations that allows to write 
down an expression for these multiloop correlators at any given genus in terms of a special 
set of Feynman diagrams. The various quantities involved depend only on the spectral curve 
of the matrix model and in particular one needs to evaluate residues of certain differentials 
at the branch points of the spectral curve. 

This algorithm and its extension to calculate higher genus terms of the matrix model free 
energy [14] represent major progress in the solution of the matrix model via loop equations 
[40-43] . This is particularly important because the orthogonal polynomial approach seems to 
fail in the multi-cut (> 2) case. A nice feature of these algorithms is that they show directly 
how the information is encoded in the spectral curve. In particular, we will be able to make 
some precise statements on the double- scaling limits of higher genus quantities simply by 
studying the double-scaling limit of the spectral curve and its various differentials. 

Given the matrix model spectral curve for an s-cut solution in the form (|2.11|) the genus 
zero 2-loop function is given by 



W(xx,x 2 ) 



1 y/cr(xi) 



2(xi-x 2 ) 2 2y/^){x x - x 2 ) 2 



(3.4) 





4(xi - x 2 )^a(xi)^a(x 2 ) 4:^a(x 1 ) x /a(x 2 ) 



where A is a symmetric polynomial given by 




(3.5) 
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where 

s-2 s-1 



° * ° ^ r dx 1 

Afo) = y)A,i4 = -zJ^fe) / ^— ^ 7 c (3.6) 



1=0 j=l 



and s is the number of cuts. The polynomials Lj(x) are related to the holomorphic 1-forms 
and defined in Appendix |X] 

Lj(x)dx f . . , 1 /0 _ N 

Wj = , / ujj = b jk , j, k = 1, . . . , s - 1 . (3.7) 



The genus zero 2-loop function for coincident arguments is 

s .. T , r/ >. cr"(a;i) cr'(xi) 2 A(x 1} Xi) 

AM (3 ' 8) 



f — 
^ Ifif-r - 



+ 



^ 16(x - o-j) 2 16cr-(x - o-j) A(x - o-j) 
The other important object is the differential 



s-l 



dS 2i -i{x 1 , x 2 ) = dS 2i {x 1 , x 2 ) = . X2 ( L% , Xl - Cj(x 2 )Lj(xi) ) dx x , 

y 'a{x l ) x 2 ^a{x 2 ) J 

(3.9) 

where i — 1, . . . , s and 

„ , x /" dx 1 . . 

CZ, a*) = / -== ? , . (3.10) 

A crucial aspect of the one-form ()3.9|) is that it is analytic in x 2 in the limit x 2 — > 024-1 or 
2 * [13] 



lim = _J_ f _J g Lj(xi) / g_ 1 ^ ) ^ . (3.11) 



The subtlety is that in the definition of (|3.1U|) . the point x 2 is taken to be outside the loop 
surrounding the j-th cut, whereas in ()3.11|) . x 2 is inside the contour. Note also that 

A(xi,x 2 ) = (^LjMCjWj (3.12) 

and in particular 

y4(xi,cri) = £j(si)cr'(crj) . (3.13) 
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The expression of W^ 9 '(xi, . . . , x p ) can be found by evaluating a series of Feynman dia- 
grams of a cubic field theory on the spectral curve [13]. To this end, define the set of 
all possible graphs with n external legs and with g loops. They can be described as follows: 
draw all rooted skeleton trees ( trees whose vertices have valence 1, 2 or 3 ), with p + 2g — 2 
edges. Draw arrows on the edges from the root towards the leaves. Then draw in all possible 
ways p—1 external legs and g inner edges with the constraint that all the vertices of the whole 
graph have valence three, namely that are always three and only three edges emanating from 
any given vertex. Each such graph will also have some symmetry factor [13]. 

Each diagram in then weighted in the following way. To each arrowed edge that is 
part of the skeleton tree going from a vertex labelled by x\ to a vertex labelled by x 2 
associate the differential dS(x\,x 2 ) (|3.9|) . To each non-arrowed edge associate a genus zero 
2-loop differential G(xx,x 2 ) = W(x\, x 2 ) dxi dx2 and to each internal vertex labelled by x% 
associate the factor (2eNy(xi)dxi)~ 1 . For any given tree T G Tp 9 \ with root x\ and leaves 
Xj, j — 2, . . . , p and with p + 2g — 2 vertices labelled by x' v , v = 1, . . . , p + 2g — 2, so that its 
inner edges are of the form v\ —> v 2 and its outer edges are of the form v — > j, we define the 
weight of the graph as follows 

P+2(/-2 2s 

W(T) = (e N y+^ II E Res ^^7^7 II dS iv {<M 

K ' v=\ i v =l ,yy VJ v inner edges v^w (3 14) 

x ] [ G 2 {x' v ,,x' w/ ) ] \ G 2 {x' v ,Xj) 

inner non-arrowed edges v'— >w' outer edges v—*j 

In order to find an expression for F g , g > 1, one should consider the same graphs relevant 
for W^ 9 '(xi) and do then the following [14]: 

(i) Eliminate the first arrowed edge of the skeleton tree. Labelling the first vertex by x\ and 
the second vertex by x 2 , this amounts to dropping the factor dS(xi, x 2 )- 



(ii) The factor (2eNy(x2)dx2) 1 has to be dropped and replaced by 

f X2 y(s)ds 

Jq ° ( J ■ (3-15) 
y(x 2 )dx 2 

Note that when evaluating the final residues at x 2 = crj, one needs to expand the above 
integral by setting q Q = o~i [14]. It is also understood that the evaluation of the residues 
starts from the outer branches and proceeds towards the root. This procedure does not 
apply for the genus one free energy whose expression has in any case been found via the loop 
equations in [44-46]. 
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We will consider the a — ► limit of each element in (|3.14j) . In this respect, it is useful 
to choose a new basis of 1-cycles {Ai, Bi}, i = 1, . . . , s — 1, which is specifically adapted to 
the factorization E — > £_ U £+. The subset of cycles with z = 1, . . . , [n/2] vanish at the 
critical point while the cycles i = [n/2] + l, . . . , s — 1 are the remaining cycles which have zero 
intersection with all the vanishing cycles. Using the results in Appendix El for the scaling of 
L, it is straightforward to argue that for a branch point bi in the critical region 

dSi(xi,x 2 ) — > dSi(xi,x 2 ) = 

U (3-16) 



where dSi is the analogous differential on S_, Eq. (J2.16j) . and for j < 

[n/2]. Conversely, the differentials dSi(xi,X2) where i labels a branch point of the spectral 
curve that remains outside of the critical region give a vanishing contribution in the double- 
scaling limit. Likewise using equations (|3.4jl . ()3.5|) . (|3.6j) and ()3.12|) we have 

G(xi, x 2 ) = W(xi, x 2 ) dx\dx 2 — >■ G(xi, x 2 ) — W(x±, x 2 ) dx±dx 2 , (3-17) 

where G(xi,x 2 ) is exactly the 2-point loop correlator on £_. 

So far we have seen that the double points of the near-critical curve do not play a role in 
taking the limit of the differentials. However, this is not the case for the final two elements 
of the Feynman rules 

ydx — ► VCa m+n/2+1 y_ dx (3.18) 

and 

I g x y( s ) ds _^ J~y-(s)d§ 

y(x)dx y^(x)dx 

To summarize: what we have found is that all the relevant quantities reduce to the analogous 
quantities on the near-critical curve in the limit a — > 0. In particular, being careful with the 
overall scaling, the genus g free energy has the limit 

F g — ► C i-s A 2-2 9 . (3.20) 

where we have emphasized that _F S (£_) depends only on £_. This is the result advertized in 
(I2.20J) and the property of universality. Similiarly, the genus g p-point loop functions have 
the limit 

W g (x 1} . . . , x p ) dxi ■ • ■ dx p — > c i-g-P/2 A 2-2g- P w g {x u . . . , x p ) dx ! • • • dx p . (3.21) 
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4 The genus one matrix model free energy 



In this section, we will consider the double-scaling limit of the genus one free energy F\ 
in more detail. This term gives information on the states that become massless at the 
singularity [7, 19]. The genus one matrix model free energy has been studied in the context 
of the Dijkgraaf-Vafa correspondence in [45,46]. In particular, the authors of [46] proposed 
an expression for a general multicut matrix model solution based on conformal field theory 
arguments by Kostov [47] and Moore [48] . This was later proved by Chekhov [44] by means 
of the matrix model loop equations. See also [45,49] for an expression of the matrix model 
genus one free energy inspired by the correspondence with the topological B model. The 
general expression is given by 

\fc=l \l<j<i<2s / V ' J ' J 

where s is the number of cuts of the matrix model solution and 

x i-i 



Nij= / —j=dx i,j = l,...,s-l (4.2) 

J A, y/<T(x) 

are periods of the holomorphic one-forms x% l^L on the reduced spectral curve. This formula 

y/<r(x) 

was derived in [44] by considering the genus one 1-point function W^ 9=1 \x) and explicitly 
inverting the relation 

A-( x ) Fl = W^ =1 \x) . (4.3) 
av 

In the previous section, we have seen that in the double-scaling limit 

W^ 9=1 \x) -> ^W {1 \x) (4.4) 

where (x) is the genus one one-point function relative to the near-critical spectral curve 

y 2 _{x) = Z m {x)B n {x) . (4.5) 
We can actually absorb the factor of A in the definition of the curve itself 

yl = tf~Z m (i)B n (x) . (4.6) 

Thus we obtain 

-4r{x)F 1 = W^ 9=1 \x) -> W^{x) = . (4.7) 

dV dV 
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We can relate the loop insertion operator ^f(^) to ^(^) as follows. Thanks to the identity 
[14,44] 

(4 ' 8) 



we find 



d(Ji , „ a / 1 \ d(Ji d<Ji 



Ax) -> + . . . )dx = a\l = a i7>(5) = "w>(z) 



This implies that 

- £(*>■ (4-9) 
Therefore, by (|4.7jl . we conclude that in the double-scaling limit 



\ 8=1 \l<j<i<" / 7 
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(4.10) 



where iVy are periods on the near-critical spectral curve S_. This is strictly correct only 
modulo the addition of a constant, but this plays no role when one considers general corre- 
lators obtained from F\ like W^'(x) in (|4.3jl . We also see that the double-scaled free energy 
depends in general on the structure of the near-critical spectral curve. In this respect, ob- 
serve that the general expression of F\ , Eq. (|4.1jl , depends on the basis of A-cycles we choose 
on the spectral curve. Upon a change of basis, which would correspond physically to an 
electric-magnetic duality transformation, F\ changes non-trivially. The expression (|4.10|) 
contains an implicit choice of basis in which the degeneration of the original spectral curve 
E into S + U £_ is made manifest [1,6]. In particular, as in Sectional we choose a basis such 
that [n/2] of the starting A-cycles shrink at the singularity and reduce to A-cycles on the 
near-critical spectral curve E_. 

In the case of the A n _i singularities studied in [1], where the near-critical curve is 

yt = B n {x) (4.11) 

we find that in the limit A( n ) — > 

71 

Fi - ~ Ya 1o S A M • ( 4 - 12 ) 
In particular, for the conifold singularity, n = 2, we retrieve the well-known result 

Fi = ~^ logA (2) . (4.13) 
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In the case of an "old matrix model" singularity, where m double zeroes collide with one 
branch point of the reduced spectral curve, <7o, and correspondingly the near-critical spectral 
curve is trivial, Eq. (j4.10|) yields 

^-—logA^o) (4.14) 
which is indeed consistent with the result given in [40]. 

The divergence of Fx and equivalently of the topological B model free energy F toPt x in the 
limit A — > indicates that there are states in the field theory that become massless at the 
singularity [7,19]. Consider type IIB string theory compactified on a Calabi-Yau space in the 
proximity of a conifold singularity. Vafa argued that ()4.13j) is consistent with the appearance 
of a single massless 4d M = 2 hypermultiplet in the low-energy theory [7]. Similarly, type 
IIB compactified on a Calabi-Yau in the proximity of an A n _i singularity yields a 4d M = 2 
theory close to an Argyres-Douglas point where mutually non-local electric and magnetic 
particles become massless [50-52] . These extra massless particles also appear in the M = 1 
theories studied via the Dijkgraaf-Vafa matrix model [1,6,38]. 

Vafa also made a proposal about a general expression for Fx 

F\ = Fx+F x = — ^ Yl lo S m * ( 4 - 15 ) 

BPS states 

where the sum is over BPS states of the M = 2 4d theory. These states can be electrically 
and magnetically charged and come from D3-branes wrapping a supersymmetric 3-cycle Ci 
in the Calabi-Yau. Their mass is given by 

m < = tJfw ■ (4 - 16) 

This is a generalization of the conifold result ()4.13|) where m 2 = |A| 2 and, as stressed in [7], 
it might not be the full answer. For A n _x singularities with n odd, the genus one expression 
(I4.12|) is indeed not matching the proposal (|4.15j) . This is probably due to the fact that 
the states becoming massless are mutually non-local. It would be interesting to understand 
better the nature of this result. 



5 Evaluation of the genus 2 free energy 



For one-cut matrix model solutions the method of orthogonal polynomials allows to evaluate 
the matrix model free energy at all genera [10]. In the case of multicut solutions, this 
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X2 




Figure 2: Diagrams (I) and {II). 



technique is not generally available. In order to evaluate higher genus terms in the free 
energy, we will resort to the recently developed algorithms that provide an exact solution 
to the matrix model loop equations [13,14] that we reviewed in Section In particular, we 
will find the expression for the genus two free energy in the case where the spectral curve 
has no double zeroes. This is relevant for the A n _i singularities considered in [1]. 

The genus two ()5.18|) and the genus one ()4.1Uj) results show that the double-scaled free 
energy depends in a complex way on the details of the near-critical spectral curve. However, 
we will see that in the simplest case, namely the conifold singularity, this dependence is 
trivial and the expressions simplify drastically. This is due to the fact that it is possible to 
choose a basis of A-cycles on the original spectral curve, before the double-scaling limit, in 
such a way that the the near-critical curve has genus zero, it is a Riemann sphere. Thus we 
recover the known result [9] 

u 2 + v 2 + y 2 + x 2 =n -> F 2 = -— (5.1) 



The explicit evaluation of F 2 involves calculating three Feynman diagrams (see FigsO 
and 01). Diagram (I) is equivalent to 

f dx 3 dx 2 dx 1 f qo S V( s )ds dS(x 3 ,x 2 ) jxr , .dS(x 2 ,xt) . 

/ = / o O O ■ ( \A O ( \ W X 3,^2 f s W{Xi,Xi) 

dxzdx. Q y(s)ds dS{x X ) 

C x 3>C x 2 2%i2in y{x 3 )dx 3 2y(x 2 ) 
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Figure 3: Diagram (III)- 

Diagram (77) is 

f dx 3 dx 2 dx x IJ y(s)ds dS(x 3 ,x 2 ) dS(x 2 ,x 1 ) . 

( n )= / . , r-j ^-^^ , v W(x 3 ,x x )W(x 2 ,x 1 ) 

J&*>c**>c*i 2m 2m 2m y\ x 3)dx 3 2y(x 2 ) 2y(x 1 ) 

dx 3 dx 2 fg S y(s)ds dS(x 3 , x 2 ) 

W(x 3 , x 2 , x 2 ) . (5.3) 



C X 3>C X 2 



2m 2m y(x 3 )dx 3 2y(x 2 ) 



Similarly, diagram (III) gives 



, dx 3 j qo y{s)ds r dx 1 dS(x 3 ,x 1 ) f dx 2 dS(x 3 ,x 2 ) 

( III ) = • , \ , / o ■ O ( \ W K X ^ X l) / o ■ O ( \ W ( X * X V 

27u y(x 3 )dx 3 Jcxi 2m 2y(x 1 ) J CX2 2m 2y(x 2 ) 



Finally, as shown in [14] 



x 3 2m y(x 3 )dx 3 



2F 2 = 2(1) + 2(11) + (III) 



(5.4) 



(5.5) 



In the following, we will only consider cases where the spectral curve has no double points, 
and we will set 

(5.6) 



y 2 = e 2 a 2s (x) 



xi n \ x 3), (Hi (EH, we find 



(I): Using the expression of W^'(x) evaluated in the appendix (jC.l|) and the differentials 



I d ^l dS ^) W(x 3 ,x 2 )wW(x 2 ) = Yl 
J c , 2 2m 2y(x 2 ) 1 3 ' 2) 2e 



1 G(x 3l a i )xf ) (x 3 ) 



16a 



1 d a" 1 

-G(x 3 , x 2 ) \ X2=ai - ^r^G(x 3 , a*) J x? Os) 



16a', dxi 



32a'/ 
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/3<7f-2<7X' v 1 d 2 ,. < Nl \ (1) , . 

+ ^ 192of G (^) + ——G^s,^)!^ - __G(x 3 ,x 2 )|„ 4 j (-3) 



+^G(x 3 ,a t ) X f\x 3 ) + (5^^3^ 2 )U 2 = CTi - |^<7(*3,*) + B^G^,.,)] ^(a*) 

(5.7) 



where we introduced 



<7(*3, x 2 ) = a ' {X3 \ + A{X ;> X2) , (5.8) 

A(x 1 ,x 2 ) = f: C ^^ , (5.9) 

C?\x 2 ) = -Y l L j (x 2 ) f p=- 1 (5.10) 

Before we illustrate how to perform the final integration, let us introduce the following 
notation 



G(x 3 ,a k ) = —==G(x 3 ,a k ) . (5.12) 



We find that 

1 if, r j„ n n\ft(s)dsJL r 1 



( ] f dx 3 J Ak Va{s)ds 

^ 2 {^J c , 3 2ni a{x 3 f/ 2 ^ 

J3<jf-2<j' i <j'>' / ,, 1 d 2 ~, a'l d ~, \ (1) 

+ [ l92a , 3 G{x 3 ^) + ——G^.x,)!^ - ^_G(*3,s 2 )|„ 4 j Xi W 

+G(x 3 , aO^^W + £ (^X? V.) + B,*JV«)) %^X, (1) (a:3) 
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+^G(x 3 ,a i )x?\x 3 ) + (^^-G(x 3 ,x 2 )\ X2=ai - ^G{x 3 , a,) + B t ^G(x 3 , a,)) tf\x s ) 
u\ \a\dx 2 2of a\ ) 

(5.13) 

The next step is to expand 

a(x 3 )^ 

in the proximity of the branch point a k itself [14]. Setting e = x — a k , we find 



where 

_3o£ 3(8< 2 - 5«') 60af - 76a' k a k a k " + lO^M'" 

Cl ' fc ~ 5< ° 2 ' k ~ 70< 2 C3 ' fc ~ 315< 3 

Finally 

(f) = 1 1 /^ c 3,fc< c 2;fc A(q- fc ,q- fc ) Ci^G^ 1 ' ^,^) -off) 



6^°)K,a fc ) + 12^4 3) K)-ar \ J_ ( 3c 3 , k a' k + c^a'j 
48 ; 16< V 4 

c li fe(6^Vfe, fffc ) + at + 18^4 2 V fc )) A^fo,^) <4V fc ) + 3<£( 2 )(a fc )' \ 
24 4 4 i 



5 fe <' \ (c 2 , k a' k c ltk A(a k ,a k ) 6A^(a k ,a k ) - a> k " + 6a> k £ k 2 \a k ) 



' 1 o\. 32< 2 ) { \ ' 1 + 2 ! 

.'S fc < \ / 3cgK , a' k ' + 3a' k C k (a k ) QA^(a k , a k ) + o' k " + Qa' k X k (a k ) + 18a' k C' k (a k ) 
' \ ~i ,9 — : r Ci fc h 



a' k 32a' k 2 J V 4 ^ 4 24 

1 f 5c 3 , k a' k 2a'l + hA{a k ,a k ) aj + 4a' k X k (a k ) + lO^^K) 
+ 32< + C2 ' fc 2 + Cl ' k 4 

QA^(a k , a k ) + aj_ + Qa' k "C k (a k ) + 2Aa' k X k (a k )' + 30a' k C'l(a k ) \ 

24 J 



192< 3 2<x£ 2 a' k 16a: 
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( ci.k<?'k A(a k ,a k ) \ 
{ 4 + 2 | 



(II): Using the expression of W(x 2 , x 2 , x 3 ), Eq. (jC.3|) . we find 

dx2 dS(x3, X2) 



C x 2 



2ni 2y(x 2 ) 



W(x 3 ,x 2 ,x 2 ) 



2s 

E 



=1 



tftPWxFW + (fe> - ^) xf'fekf (,3) 



+ 



AW{<Ti,<Ti) , A(a,,a,) 



8a' 



+ 



Ifi n-'f 



16a' 2 



(5.14) 



xf^s) 2 



16(cr j - Oi 



a' 2 + a> 2 



16a>a> 



16(ct-< — o-j)^- 



(5.15) 



Similarly for (//) we find 

2s , r „(3) 

"l6 



1 6£V fc 16 

fe=l K 



C 1)fc A((T fe , CT fc ) + A (1 ' 0) (cr fc , (T fc ) 



/ A(q fc ,q fc ) < \ / (2) 



/ AW(a k ,(T k ) A(a k , a k f _ ^ y- 1 _ <r k A{a k , a k ) \ f 2A(a k , a k ) 



C7 



+E 2 



' 2 +<rf\ , A(a k ,a j f , Afa, ^ - ^)\ / 1 %,a 



+ 



16cr^.crj 16(o- fc - a^a'f.a'j J \a k -(r 



(5.16) 



(III): From Eqs.flED and (jSZQ) , we find 



fc=l K 



c 3 , fc + 2ci, fc 4 2) (a fc ) + 24 z;/ ((7 fc ) 
256 



(2)/ 
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Cl,fc 

128 



E 



+ /f (<*)) + 



128 



E 



((7 fc - a^) 3 



+ £ (2) 'K) + 



8 



+^(c 2 , fc + 4 2) K) + 



Cl,k 

8 



(Tfc - (To 



+ 



i4(<Tj,<T fc ) 



+ 



8 V (o-fc - ^ 



+ 



^4 (o-fc, o-j' 



(5.17) 

Therefore, the final expression for the genus two free energy is 



— / 157af 491<af 35af" 35af A(a,, a,) lltrfAfo, a*) 49afA(a i , a,) 2 



£ 2 4f V 15360a' 4 46080a' 3 3072a' 2 768a' 4 576a' 3 640a 

5A(a,,a t ) 3 _ 37af£( 2 )(a t ) A(a t , a,)A (2) (a,) 13£( 2 )'(a t ) 54 3) (a t ) _ 11<£^) 
96af " 2560af * " 24af ~ 1536a,' 384a^ " 768a,' 2 

A(a ti a^(a t ) 5£?fo) _ 47q^ 1 )(q t , a,) 5A(a,, a,)^ 1 )(a t , aj ^ 2 )(a t! a t ) 
32af 768a^ 7680af 384af 768af 

89a l "A( 1 '°)(a l , a,) a^ 1 - )^, a,) + 7A(a„ gQ^fo, a t ) + ^(a,, a t ) + 5A( 2 '°)(a l? a, 



3840a' 3 160a' 2 96a' 3 384a' 2 768a' 2 



i 

2s 2s 



^ ZS AS / / 

z 2 U j£ V 768(a, - a,)^ + 48(a, - a,fa] + 48(a, - a,) 3 af 
a!' a" aj afa" 



384(a, - a,) 2 a 4 2 + 10(a, - a,) 2 a^. + 1536(a, - a,) 2 a^ + 384(a, - a,)af a^ 

^fo,Qj) _ ^fa, g») aJA(a^) 
128(a, - a^af 3(a, - arfa'a' 128(a, - a>^ 

Afaaj) Ma^oj) A{a i ,a j ) a'l^A{(T i ,a j ) qVg'!A{a i ,aj) 



24(a,-a,) 2 af 48(a, - a,) 2 af 48^ - a,) 2 a^ 512afaf 1536afa^ 



a'JAiai^^Aia^cTj) a^A{a i ,a i )A{a i ,a j ) Afaaj) 2 A(a h a A 2 A^^a^) 



192afaf 384afa^- 48 (a* - a^af ' 24(a^ - a^)af aj ' 48afa;. k 
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384(a l -a,) 2 ^ 96^ - a^af a] 32(^ - <7>f^ 



384<rf^ 48(xfaf 96af(xj 



512af 192(rf 1536^' 1536^'af 384^af 

Let us consider the limit 

y 2 = e 2 a 2s (x) — > e 2 (x s - a s ) , a — > 
where s of the branch points come together. In the double-scaling limit 

a ^ , e ^ oo , A = ea s/2+1 = cnst (5.20) 

we find 

F 2 — ^g^, E_:id = ai'-a' (5.21) 

as explained in Sectional However, the final result will not simplify in general. In fact, it 
depends on the details of the near-critical spectral curve, which has genus [(s — l)/2]. An 
exception to this is given by the case of the conifold singularity, where the original spectral 
curve becomes in the limit 

y 2 ss e 2 {x - a)(x - b) , a, b -> , (5.22) 

which is essentially the spectral curve associated to a Gaussian matrix model [6]. As in 
Sectional it is convenient to choose one of the A-cycles of the original spectral curve to be a 
loop encircling the cut going from branch point a to branch point b. This cycle will reduce 
to an A-cycle on the near-critical spectral curve 

y 2 _ = a(£) = (x-a)(x-b) (5.23) 

where, as before, the tilded quantities are finite in the limit a, b — > 0. The above near-critical 
curve is actually a Riemann sphere. In particular, one can check by evaluating the residues 
of all the integrands at infinity that all periods of the form 

>iV 1 (5.24) 



'a(x) [x - &i ) n 



are zero. The expression of F 2 simplifies dramatically 

4 



15e 2 (a - bY 240S 2 
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(5.25) 



where 

S = fydx-+?£^-~A, (5.26) 

and in terms of S the genus zero free energy is given by 

1 ap i 

F ° W 2 S ^i W 2 52l0g5 ' (5 - 27) 

Thus, ()5.25|) indeed matches the expected result for the genus two free energy at a conifold 
singularity [9], which is equivalent to the c = 1 non-critical bosonic string [20]. This par- 
ticular singularity is obtained from a 2-cut solution with a cubic superpotential in the limit 
where the two cuts touch each other. The fact that this limit should be equivalent to the 
c = 1 non-critical string was also observed in [23]. 



6 Conclusion 



The class of matrix model DSLs that are associated to the large N field theory DSLs intro- 
duced in [1] define a class of c < 1 non-critical bosonic strings [6]. They fall into different 
universality classes from the ones usually considered in the old matrix model. We argued 
that these non-critical bosonic strings are related to the topological twist of non-critical 
superstring backgrounds of the form SL(2)/U(1) x LG(X n ) that are dual to the large N 
double-scaled field theory and the associated four- dimensional double-scaled LST at the 
corresponding A n _i singularity. To study the matrix models, and the relevant multicut solu- 
tions, we used the techniques of Chekhov and Eynard based on loop equations. These allow 
to show in general that the scaling of the higher genus terms in the perturbative expansion 
of the matrix model free energy matches precisely the scaling of the topological B model free 
energy in the vicinity of the Calabi-Yau singularity, which is consistent with the Dijkgraaf- 
Vafa correspondence. We also evaluated the genus one and two terms explicitly for the A n _i 
singularities, recovering the conifold result in the n = 2 case. These techniques allow to 
study multicut solutions where the "old matrix model" tools are not generally available, 
but further work would be needed to find the exact expression of the perturbative matrix 
model free energy at all orders for the A n _i singularities with n > 2. In particular, it would 
be interesting to see if this perturbative series needs a non-perturbative completion like in 
the conifold case. Such completion should correspond to D-brane effects on the non-critical 
string side as in [53]. It would also be interesting to perform the topological twist of the 
SL(2)/U(1) x LG model and determine the non-critical bosonic string explicitly. 

Acknowledgments. I would like to thank Nick Dorey, Tim Hollowood, Luis Miramontes, 
Sameer Murthy and Asad Naqvi for various discussions and comments. 
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Appendix A: Some double-scaling formulae 



In this appendix, we consider the double-scaling limit of various quantities defined on the 
curve £ (j2.11j) . This is most conveniently done in the basis {Ai, Bi} of 1-cycles described in 
Section EJ In particular, for i < [n/2] these are cycles on the near-critical curve £_ in the 
double-scaling limit. 

The key quantities that we will need are the periods 

Mij = (t> — dx , Nij = (t> dx . (A.l) 

First of all, let us focus on iVy where j < [n/2], but i arbitrary. By a simple scaling argument, 
as a — > 0, 



N l3 = I U) 4== dx — a-/ 2 / U) 4= dx = a-/ 2 /f >(&,) , (A.2) 



for some function ' of the branch points of S_. Here, are the two branch points 
enclosed by the cycle Aj. A similar argument shows that scales in the same way: 

M^-^a^f^Kh) . (A.3) 

So both and M^, for i, j, < [n/2], diverge in the limit a — > 0. On the contrary, by using 
a similar argument, it is not difficult to see that, for j > [n/2], and are analytic as 
a — > since the integrals are over non-vanishing cycles. 

In summary, in the limit a — > 0, the matrices N and M will have the following block 
structure 

Nfl ' * ( Mf I ' (A ' 4) 



where by — or + we denote indices in the ranges {1, . . . , [n/2]} and {[n/2] + 1, . . . , s — 1} 

respectively. In (jA.4|) . iV and M are divergent while the remaining quantities are finite 

as a — > 0. 

We also need the inverse L = AT -1 . In the text, we use the polynomials Lj(x) = 
Ylk=i LjkX k ~ l ", which enter the expression of the holomorphic 1-forms associated to our basis 
of 1-cycles, 

Uj = Sij . (A.5) 
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These 1-forms are equal to 

Uj (x) = ^= dx = EL ^ fc-1 dx , I cvj (x) = 8 tj (A.6) 
where i,j — l,...,s — l. From the behaviour of N in the limit a — > 0, we have 

L = iV- 1 — ( Y ^ -i j , AT = -AC ! ivi°> (a^) _1 . (A.7) 

Since A is singular we see that L is block diagonal in the limit a — * 0. This is just an 

expression of the fact that the curve factorizes S — > XL U S + as a — > 0. In this limit, using 
the scaling of elements of L jfc , we find, for j < [n/2], 

cjj — >■ ==i = cjj . (A. 8) 



the holomorphic 1-forms of E_. While for j > [n/2], 



Uj — ^ fc>[ " /2]( 7_^ rfx , (A.9) 



are the holomorphic 1-forms of E+. 



Appendix B: The explicit expression of Xl (p) 



Using the formalism developed in [13] and [14] to solve the matrix model loop equations, 
one can easily find the expression of the differentials Xi^ (p) defined by 

(k-2W (p)) X < t\p) = ¥ ^, (B.l) 

where <7j is a branch point of the matrix model spectral curve. These 1-differentials appear 
quite naturally in the expression of higher loop correlators in the matrix model and in the 
integration steps leading to F 2 . We have 

v {n) (p) = Res ( dSiM - ^ (B 2) 
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Given the expression of dSi(p, q), we can easily perform a Taylor expansion in q around the 
branch point <7j and find the residue. We will mainly consider the case where the matrix 
model spectral curve has no double points, setting 

y 2 = e 2 a 2s (x) . (B.3) 

In this case 



d Sjjp, g) = dSj{p,q) = 1 / 1 _ L , v f dx 1 
V(<l) ~ ~ e^P) \p-q jr[ iKP) J Aj ^{x) 



(x - q) 



dp 



Note also that the expression in brackets is analytic in q. Then, for instance, we find that 

xf\p) = — = i — - > ; i - , i dp 



2e^P) ip-Vi |J Ljip) Sa 3 



X) [X-CTi) 



' ( 1 +£ i (p))dp=-^=(— + ^^)dp (B.4) 



2ey/a (p) \P ~~ a « / 2e^/a{p) \P - d cr[ 

Xt\p) = 7= ( 7—^ -YHP) I ^^7—^ 1 * (B.5) 

2e^/aipj\(p-Vi) 2 jr[ ~ Ja 3 y/^xj(x-ai) 2 J 

and in general 

1 1 d n ^ ( 1 s ^ f dx 1 \ 

The above expressions can be generalized to the case where the spectral curve is of the form 

y 2 = M(x) 2 a(x) , (B.7) 



(n), v 1 1 d^ 1 1 / 1 ^ T , , f dx 1 V 

X * ^ " 2v^) (n - 1)! dq^ M(q) [p - q £ ^ ^ ^ (* - ,) ) ' 

(B.8) 

Appendix C: Evaluation of W <yl \p) and W(p,p, q) 
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In this section, we are going to evaluate two loop-functions whose expression is needed 
later for F 2 , the genus one one-loop function W^(p) and the genus zero three-loop function 
W(p,p, q). Let us start from W^(p). Using the diagrammatic rules of [13] we find 



2s 

> 2 ) = Res Xl ^ a% 
i=i 



2s 



1=1 



dSi{X2, X\ t 
2y(xi) 

1 



Bi 



dSi(x 2: xi) 
2y(x 1 ) \l6(x 1 - (Ti) 2 ' x 1 -a i 



i=l 



(C.l) 



where 
B 



o"{<Ji 



<j"{<Ji) + A(ai,ai 



(C.2) 



16<r'(<7i) 4a'(a t ) 

This is exactly the expression given for instance in [44], once we use the identity 

Then, let us evaluate the genus zero 3-loop function with two coincident arguments W(x2, x 2 , x 3 ). 
We find 

' dSi(x 2 ,xi) 



W(x 2 ,x 2 ,x 3 ) = ^2Res xl ^ ai ( 
i=i ^ 

2s 



i=i 



dSi(x 2 ,Xi) 



2y(xi 



W(xi,X2)W(xi,X3) 
A(x 2 , (7, 



+ 



2y(xi) U^-^V^N 4v/^2" 



4(x 3 - <Ti)y/a(x^ 4^<r(x 3 ) / <7'(^)(^i - <T< 



2s 

E 



A(x 2 ,<Ti) 



v 4(x 2 - (Ti)y/a(x 2 ) 4 v / ( j(x 2 ) / \4(x 3 - ai)y/a(x 3 ) 4 v /cr(x 3 ) / ^VO 

-\r(.r 2 )M- J '(.r:0 ■ (C.3) 



2s 2 // \ 

4 



i=i 
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